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The boundary-shape formalism of Morse 
and Ingard is applied to the acoustic 
modes of a deformed spherical resonator 
(quasisphere) with rigid boundaries. 
For boundary shapes described by 
r = a[\- eT{Q, (^)], where e is a small 
scale parameter and :^ is a function of 
order unity, the frequency perturbation is 
calculated to order e^. The formal results 
apply to acoustic modes whose angular 
dependence is designated by the indices 
/and m. Specific examples are worked 
out for the radial (/= 0) and triplet (/= 1) 
modes, for prolate and oblate spheroids, 
and for triaxial ellipsoids. The exact 
eigenvalues for the spheroids, and 
eigenvalue determined with finite-element 
calculations, are shown to agree with 
perturbation theory through terms of order 
E^. This work is an extension of the 
author's previous papers on the acoustic 



eigenfrequencies of deformed spherical 
resonators, which were limited to the 
second-order perturbation for radial modes 
[J. Acoust. Soc. Am. 71, 1109-1113 (1982)] 
and the first order-perturbation for 
arbitrary modes [J. Acoust. Soc. Am. 79, 
278-285(1986)]. 
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1. Introduction 

Spherical acoustic resonators [1] have been success- 
fully applied to measurement of the universal gas 
constant [2] and to gas thermometry [3-7]. The radial 
(0, fi) acoustic modes are well-suited for high-accuracy 
work because they are non-degenerate, well separated, 
insensitive to visco-thermal boundary effects, and only 
weakly dependent on the details of small shape imper- 
fections [1,8,9]. Nearly-spherical resonators, now 
referred to as quasispherical resonators [10,11], have 
been designed to facilitate the measurement of the 
electromagnetic resonances for determination of the 
mean radius of the quasisphere. The speed of sound can 
be determined from the combination of measured 
acoustic and electromagnetic resonance frequencies 



[12]. The acoustic and electromagnetic eigenvalues for 
quasispherical resonators must be evaluated using 
approximation methods. In this paper, boundary-shape 
perturbation theory [13] is used to calculate the 
acoustic eigenvalues. 

An idealized rigid spherical resonator of radius a has 
acoustic modes with the acoustic pressure proportional 
to the eigenfunctions 



<!>.„„ =7. (^.„'-)^rJ^,<A), 



(1) 



where 7/ ((^ is a spherical Bessel function and Y^^ is 
a linear combination of spherical harmonics 



^ ^mii^l^ 



(2) 
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with coefficients a^^ chosen to make the Y^^ real. The 
eigenvalues for the perfect spherical geometry are 
Kn = ^in /«, where (^^„ is the nth root ofy^C-^J = 0. 

The bounding surface of a quasispherical resonator 
has the form 



:a[l-£j^(0,0)], 



(3) 



where £ is a scale parameter satisfying < £ « 1 and 
j^ is a smooth, non-negative function of and (j). The 
£ = eigenvalues of a quasisphere differ from those of 
a perfect sphere of the same volume by a fraction of 
order £^ or higher [1,8]. The non-radial acoustic modes 
of a perfect sphere occur in multiplets of degeneracy 
2^ + 1 . Typically this degeneracy is split to order e in 
a quasisphere, but the mean eigenvalue of any multiplet 
differs from the corresponding eigenvalue of a perfect 
sphere of the same volume by an amount of order £ ^ or 
higher [9] . The same is true of the electric and magnet- 
ic modes of a quasisphere with perfectly conducting 
walls [12]. 

In principle, a highly accurate measurement of the 
speed of sound in a gas can be made by measuring the 
acoustic and electromagnetic resonance frequencies of 
the same quasisphere. Geometric contributions to the 
error will then be of order e ^. If the shape is known, 
higher accuracy can be obtained if the theoretical coef- 
ficients of the £^ perturbation terms can be calculated. 
This has already been achieved for the radial acoustic 
modes [8]. A corresponding theory for the electromag- 
netic modes, a subject of current research by the author, 
is much more complex, and is closely related to the 
theory of the second order shape perturbation theory for 
the non-radial acoustic modes, as developed in this 
paper. The results derived here will be useful in exper- 
imental studies of quasispheres which will compare the 
effects of shape on the acoustic and electromagnetic 
spectra. Also, the results enable the use of the lowest- 
frequency acoustic modes, the 11m triplet, to be used 
for high-accuracy work. 

2. Formalism 

Morse and Feshbach [13] (MF hereafter) present a 
formalism for calculating the eigenfrequencies of an 
acoustic cavity resonator C enclosed within an unper- 
turbed cavity Cq. Figure 1 illustrates the geometry; the 
surfaces enclosing C and Cq are designated S and Sq, 
respectively and the region between C and Cq is desig- 
nated C. The unperturbed cavity has a set of eigenfunc- 
tions 0;v and eigenvalues k^ , satisfying the Helmholtz 
equation 



in Co and the Neumann boundary condition 
3^. 



dn 



" -nV^^=0 



(4) 



(5) 



on Sq. (For brevity, the subscript A'^ in these equations, 
and other subscripts in capital letters, are used to repre- 
sent sets of lower-case mode indices.) 




Fig. 1. Perturbed cavity C (boundary S) within unperturbed cavity 
Co (boundary Sq). The region between C and Q is designated C^. 

The perturbed problem is defined by a surface S 
enclosed within Sq, and enclosing a cavity C. The per- 
turbed problem satisfies 

(V'+A:')O = 



(6) 



in C and the Neumann boundary condition 



dn 



= 



(7) 



on S. A second-order expression for the perturbed 
eigenvalue Ic^ is [MF Eq. (9.2.53)] 

[A,,-{e-kl)^f,,][A,,-{l^-l^,)^f^^] 



where 



and 



j^,,{e-kiyA,^ 






(8) 



(9) 



Ko=^^QP=Jc^P^QdV. (10) 

A more useful computational form for the integrals 
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(9) can be obtained by applying the divergence theorem 
to OpVOg ± OgVOp in C to obtain 

(11) 
and 

^qp-Aq = (k'^-k',)l^^,%dV, (12) 

A negative sign occurs in these expressions because 
the outward normal from C on 5* is -n. For Eq. (12), 
use was made of the the orthogonality of the unper- 
turbed fiinctions in Q for p^q: 



J^O^OgC/r + J^O^OgC/F=0, Q^P. (13) 

An expression for the diagonal terms follows 
directly from Eq. (11): 

Ap=l[K\'^pf-\^'^pf]dV. (14) 



This is an integral over the region between S and 5*0 
of a quantity proportional to the difference between the 
potential energy and the kinetic energy The correspon- 
ding term for the perturbation of the electromagnetic 
modes has the same form [12]. An expression for the 
off-diagonal terms can be obtained from the sum of 
Eqs. (ll)and(12): 

A,,=l[kl0,%-V^,-V^,]dV, (15) 



3. Deformed Spherical Resonator 

Consider a deformed spherical resonator with a 
boundary surface S defined by Eq. (3). 

When applied to a quasisphere with S defined by 
Eq. (3), the volume C, and accordingly the integrals in 
Eqs. (14) and (15) are of order e. Equation (8) can then 
be solved iteratively to obtain 









(16) 



u„ 



^N M^N^kl(kl-kl) 



-0{e^). 



To evaluate Eq. (16) to order £^, the numerator of the 
first term on the right must be calculated to 0{£^), the 
denominator to 0(e), the coefficients A^^ in the sum 
term to 0{£). The normalization constants in the 
denominator of the sum term need only be calculated to 
0(1). 

The acoustic modes of a perfect spherical resonator 
occur in multiplets with (21 + l)-fold degeneracy. Only 
the f = radial modes are nondegenerate. When calcu- 
lating the perturbation series for nonradial modes, the 
coefficients a^^ in Eq. (2) should be chosen to make the 
coefficients A 



PQ 



zero for the modes with kp = kq. 



Equation (12) shows that A^^ = App^ exactly for degen- 
erate pairs, so the proper choice of coefficients a^^ can 
be obtained by diagonalizing the submatrix [Aj^] link- 
ing the multiplet terms. 

More precisely, the off-diagonal terms of this sub- 
matrix must be of order £ ^. Consider the application of 
Eq. (8) to the multiplet components. The first-order 
perturbation shift of each component is of order £. 
In an iterative solution of Eq. (8) the order of the 
terms in the denominator of the sum terms would be 
Npp = 0(1), I^-k^= 0(£% and App = 0(£), The numer- 
ator is the square of Aj^p -(I^-k^) Afpj^, Both 1^-k^ 
and A/Jvp are of order £, so if ^ivp = 0{£^) the numerator 
will be of order £^, and the entire sum term of order £^. 

To get the first term in Eq. (16) to 0{£^), the numer- 
ator must be calculated to 0{£^), and the denominator 
to 0{£). The thickness of the integration volume C in 
Eq. (14) is of order £, so the integrand of the numera- 
tor is needed to 0{£). Within C'the radial derivative of 
the spherical Bessel fimction is of order £ so the fimc- 
tion itself satisfies 



Ji{Knr)=Ji{^M)+0(£'y 



(17) 



Integrals Npq with Pi^Qdo not appear in Eq. (16), 
only normalization integrals for which the repeated 
indices are superfluous. The notation can hence be 
simplified by using an ordinary math font for M and 
a set of lower-case indices to designate the mode. In 
the new notation, the normalization integral in the 
denominator is 

N,^=\da[;[j,{k,„r)i\tjr^dr, (18) 

which can be evaluated as the difference between an 
integral from r = to a and an integral from r^ to a to 
obtain, for inm 7^ 1 0, 
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^inm ~ ^ Jt \}Stn )\ 






(19) 

For the special case nlm = 010, the eigenfunction is 
Joi^oif) ^ 1 and the eigenvalue is <^oi ^ 0- The normal- 
ization constant is 



A^„, 



0(8), 



(20) 



which differs from Eq. (19) by a factor of 2/3. 

The function JT may itself depend on the scaling 
parameter e; it is usefiil to make this explicit: 



(21) 





T-- 


= ^,+£^+0(£^). 


The first term 


in Eq. (1 6) then has the form 


A 

inminm 







^l-i{i^\)-lE^l\T,\X,J^ da 
+0(£'). (22) 

(Note that, the operator rV appearing in this expression 
involves only angular derivatives.) The coefficients 
(15) in the perturbation series are only needed to 0{e), 
so the integrand is only needed to 0(1); only the lead- 
ing order of Eq. (18) is needed in the denominator. The 
sum in Eq. (16) simplifies to 



Inmtnm' \ 



'n'm'^lnm ^ inm^ I'n'm'^ln y^ln K'n' ) 



tnm 



(23) 



^*^ \^iml'm'\ 






' t,r.^n,n^i-i{i+\) [0,-i'i('+i)]{a-&„') 



where 



The sum over n'in Eq. (23) is 

(25) 



°° . 1 B' 



where the prime on the summation symbol indicates 



the omission of the terms with i'v= in. The sums are 
evaluated analytically in the Appendix using the 
technique of Ref [8]. The results for -£ V are 

'^''' ^^^'" ) forf^i, 

^^JtiLY (26) 

^.t -3^(^ + forr = ^, 



and, for f = 0, 



1 I >o(4) for -£7^0, 

2a 2U0L' (27) 

-l/(4^o'«)' for^ = 0. 



The full sum in (16) is thus 

\A I' 

X^ \ inmtnm' \ 



fn'm'^Mm ^ Inm^ tnwKn {J^ln Kn' ) 



(28) 



Ae^ 



^i-K^+^)i 






2e ^Is^"^ 



^i[^i-i{i^^)Y 



where the last term is 1/3 of the contribution from the 
010-mode, which has a special normalization; 2/3 of the 
contribution of this term is included in the sum term. 

3.1 Reference Eigenvalues 

In order to separate out the effects of shape from the 
effects of volume, the perturbed eigenvalues J^ will be 
compared with the eigenvalues (^/„ J^ of a reference 
sphere of the same volume V as the perturbed sphere. 
The fractional difference equals 



y-inm ) 



£2 

■sin 



(29) 



where {ka^^l^\, is the product of {a'laf and the sum 
of 1 and the series on the right side of Eq. (16). 
The ratio of the volume V = 4K{a ')V3 to the volume 
Fq = 47raV3 of the unperturbed sphere is 



7 hilo-^)'"" 



4;r 



(30) 



= 1-3s{T,)+3s^[t^^ -T^)+0(s'). 
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where the triangular brackets indicate an average over 
solid angle. The ratio of the squared radii is 

/\2 



ffl 



■ \-2e{K)+e^ 



(TS+llT^-T, 



(31) 



The desired fractional difference is 

4£^ 



{ka)'-^i_A, 



<? 



k N 



S-^(^+l) 



Z|d(«) Pc 



28 '\B 



2 d(«) 



£mOO\ 



A 



-2£^n, 



k N 



-2£(T,)+£ 



0{e^), 



t-(^o) 



+2[T')-2{T, 



^>] 



(32) 



pressions (34) with \m'\<t and | m |< ^. Accordingly, 
it is possible to identify the terms that can possibly con- 
tribute to non- varnishing values of 5^^^^'^' , by applying 
the following rules: 

1. First look at the non- vanishing Cx^ in Eq. (33). 

2. For the unperturbed mode index i, look at each 
expansion-coefficient index X and find the 
values of/ satisfying the conditions of Eqs. (36) 
and (38). 

3. Note that the coefficients a^^^ in Eq. (2) are 
often non-zero only for //' = ±m. Consider the 
possible terms; then for each unperturbed mode 
index m, and expansion coefficient index // find 
the value of m' satisfying Eq. (37). 



where the term coupling to the 01 -mode has been made 
exphcit. 

3.2 Series Evaluation 

Identification of the contributions to the coefficients 
^tmi'm' is facilitated by expressing the shape as 






(33) 



Equation (24) then involves linear combinations of 
terms of the forms 



{rm'\Xfl\im)^JY:,^,Y,J,JQ. (34) 



and 



r 



^Y;^ dY,„ 1 dY^ 8K ^ 



>^ "■'tm 



de de sin'e 30 30 



dQ. 



1 ^' 



(35) 



where Eq. (35) was obtained using the technique 
described in the Appendix of Ref [9]. Alternatively 
Eq. (35) can be derived using the raising and lowering 
angular momentum operators (see, e.g. Ref [14]). The 
bracket expressions (34) vanish unless the following 
conditions are satisfied: 

1. \r-i\<x<r+i, (36) 

2. m' = fi + m, (37) 

3. / + A + f must be even. (38) 

It is clear from Eqs. (2), (34), and (35) that B^^].^. can 
be expressed as a linear combination of bracket ex- 



4. Once the possible non-vanishing coefficients 

^tmtm' ^^^ identified, computation of the values of 

the coefficients can be carried out using symbol- 
ic algebra software. 

4. Examples 

The second order perturbations of the i = radial 
modes and the three-fold degenerate i = 1 modes are 
worked out in this section for prolate and oblate ellip- 
soids, and for triaxial ellipsoids. 

The eigenfunctions of the unperturbed i = modes are 



^OnO=Jo{Kn^)^0^ 



(39) 



The appropriate i = 1 unperturbed eigenfunctions 
for any quasi-spherical resonator that has its major axes 
aligned with the x, j, and z directions are the product of 
jiikind and 



Y JLl^^ J_sinacos0=j— - 
'' V2 V4;r V4;r r 



^10 =^n = J COS0 = J 



(40) 



Y. .= — ''' =J — sm0cos<^=J — — , 
' V2/ V4;r ^ V4;r r 

for which the submatrix with components ^i^^j^^' is 
diagonal in mm \ 

4.1 Prolate Spheroid 

For a spheroid of semi-major axis a and semi-minor 
axes b = al{\ + £), with £ > 0, the radial coordinate is 
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with 



^l + (26;+6;')sm' 
J^ = s'm^ O+ei -s'm^ --sin' oYo(£^), 



(42) 



for which 



(^o)=f (^0^4, (^)=-^. (43) 



The shape function J^q is an exact linear combination of 
Yqq and 720- Accordingly, for the radial modes, the con- 
tributions to ^oorm' ^^^ limited to ^'=0 and f=2. For the 

i = 1 modes, the contributions to B^'^^^^^ are limited to 
r=l andf = 3. 

4.1.1 ^=0 Modes 

The non-vanishing coefficients are 



d(«) _^;:2 D(n) _ 2V^ £2 

^0000 T So«' ^0020 ^ ^ *=0« 



(44) 



Evaluation of Eqs. (79) and (78) yields 

^ s =1 

' *^0rt2 ^ ' 



4£' 



(45) 



where the latter was obtained using the condition 
Jo i^on) = -j\ i^on) = and recurrence relations for the 
spherical Bessel functions. Evaluation of the perturba- 
tion series (32) yields 



(kci'f-^i_8^le'^ 



Son 



135 



-0(e^), 



(46) 



in agreement with Eq. (30) of Ref [8]. 

4.1.2 ^=1 Modes 

Equation (22) is 
4 

4(g.^„-4) 2(94t-1264>440) 

A 

^ln±l,ln±l 



(47) 



0{e^) 



4(gf„-3) ^ 2(6^: -2594^2 70) 



e'+ 



0{e^). 



The non- vanishing coefficients are 

^10,10 ^ V?lrt ^/' 

5(«) ^_2V2T. , _ . 

^10,30 T^ V?lrt ^/' 



(48) 



^1,±1,3,±1 



2VI4 

35 



(^«-4> 



From Eqs. (74) and (72), the required sums are 



5i„i = 



4(C-2T 






10(4^„-4)' 



(49) 



the latter following from recurrence relations for 
the spherical Bessel functions and the condition 
7/ (4„) = 0. Substitution of Eqs. (43) and (47)-(49) into 
Eq. (32) yields 

^u ~ 15(4^„-2) 

4(54^^ +3734t +495^ -6924^,^„ +2980 ) ^ 



7875 (.^,^„ -2)' 
+0(e'), m = 0, 



(50) 



and 



^.^. "15(^,^2) 

2(72^f„ -961^f„ +3285^ -3282^,^„ -2560 ) ^ 
+— ^^ ^ ^£ 



7875(.^,^„-2) 
+0(€^), m = ±]. 



(51) 



The scalar Helmholtz equation separates in 
spheroidal coordinates, so the acoustical eigenvalues 
can be determined by direct numerical calculations 
[15]. The eigenvalues calculated numerically for a 
series of values of e are compared with Eqs. (50) and 
(51) in Figs. 2 and 3. 
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Fig. 2. Comparison of perturbation series (lines) for prolate spheroid 
with exact numerical solutions (symbols). 
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The average eigenvalue perturbation for the 
In-triplet is, from Eqs. (50) and (51), 



'(ka'f-^t 



£2 
31« 



(52) 



S{3C -14^L +90^;„ -243^ +10) ^^ ^ 

1125(^.^„-2)^ '^ ^' 



which has no linear term, consistent with the general 
results derived in Ref [9]. 

4.2 Oblate Spheroid 

For an oblate spheroid of semi-major axis a and 
semi-minor axes b = a/(\ +£),£> 0, the radial coordi- 
nate is 



^l + (2£+£')cos' 



= a(\-eJ^), (53) 



with 



T = cos' e +£| -cos' — cos' \+0(e^ ), (54) 



0.001 0.003 



0.01 



0.03 



for which 

(^0) = -' {^o') = -> i^^)=—^- (55) 

The perturbation calculations for the ^ = and i = 1 
modes parallel those for the prolate spheroid and will 
not be reproduced in detail here. The final expression 
for the eigenvalue perturbations for the radial modes is 
exactly the same as the result for the prolate spheroid 
[Eq. (46)]. For the £ = 1 modes the fractional perturba- 
tions are 

^ 4(5< -6IIC +3645g;„ -6924C -1220 ) ^^ ^^^^ 

7875 (^,'„ -2)' 
+0(£'), m = 0, 



Fig. 3. The absolute difference between the exact numerical eigen- 
values k\^^^ and the predictions ^^p^rtOf Eqs. (50) and (51) plotted as 
a fraction of ^^i^ = (^Ja'f', as functions of e. The dashed line, intend- 
ed as a guide to the eye, is proportional to e^. The plots show that the 
differences are approximately proportional to e\ The numerical 
values exceed the perturbation values for the 1 1 1 mode over the full 
displayed range, and for the 140 mode for e< 0.025; for all other 
cases the difference is negative. 



{ka'f-^L 4(^.>1) 

^ 2(12^1 +ml +135^ -3282^.^„ +1640) ^ 

7875(^/„-2j 
+0(e'), m = ±\. 



(57) 
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Exact solutions for the oblate spheroid [16] were calcu- 
lated and compared with Eqs. (56) and (57). The agree- 
ment, like the corresponding agreement for the prolate 
spheroid, is very good. The plots resemble Figs. 2 
and 3. 

The mode average is exactly the same as Eq. (52) for 
the prolate spheroid. 

4.3 Triaxial Ellipsoid 

The surface of the triaxial ellipsoid defined by 



- + yV- 



(58) 



(l + £,) (1+8,) (1 + 8,) (1+8,) 

can be expressed in the form of Eq. (3) with 

8J^Q =8^sm^ +82(00$^ +sin^ Osin^ <^) (59) 

and 



8^J^,=8,\-cos^e-l Isin'e 



+^i^2(~sin^ <^-3cos^0cos^<^')sin^ 

(5 3 

+^2 -1+ — sin^0cos^<^-— sin"^ 0cos 



(60) 



The shape T^ is an exact linear combination of Fqo, 
720, ai^d 72 + 2- The non- vanishing values of ^j^r,„' are 
accordingly limited to the same values of t as for the 
spheroids. 

4.3.1 ^ = Modes 

The non-vanishing coefficients are 



*^ ^0000 ^ ^n ^*^ 1 ^ *^ 2 / 5 

Evaluation of the perturbation series (32) yields 



(61) 



Son 



-^{el-e,e,+el)+0{e'). (62) 



The correctness of this result was checked by calcu- 
lating the radial-mode eigenvalues of a triaxial ellipsoid 
using a fmite-element method. The parameters 8^ and 82 
were varied, with the ratio held constant at £1/^2 = 2. 
Figure 4 shows that the difference between the finite- 
element results and Eq. (62) is cubic in 8^. 



CNO 






I -4 







(ei/0.01)^ 

Fig. 4. The differences between the values ofk^^^^ determined with 
the finite-element method for the 02 radial mode and the predictions 
k\^^ of Eq. (62) for an ellipsoid with £2 = ejl. The line proportion- 
al to £1 was fit to the plotted points. 

4.3.2 i =1 Modes 

The non- vanishing coefficients B^'^^.^. are, with ^ = ^^^ 
for brevity, 

5ro!.,o=|[(r-4>,+(2£^-3>,], 
<>,,=|[(2r-3>,+(r-4>,], 

Bl%,,_,=-^(^'-4)(4e,-e,), (63) 
<!3,o=-^(r-4)(2£,-£,). 



<k.=-^0'-4)(4£,-3£,), 

^1.0.3,2 /^^ ^^2 ' 

<k3=-^r-4K. 

The fractional perturbations for the Inm modes have 
the form 






15(r-2) 



Pnn,£l+Pl2n,£l£2+P22n,£l 
7875 (^' -2 1 



(64) 
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1 ~ P2,Q ~ P2,-l ~ ^5 

= -2, 



(65) 



where 

A,i = / 

Pi,o = I 

p,,^, = 2(72^' -961^' +3285^' -3282^' -2560) , 

p^^^^ = -8(27^' -76^' +1035^' -3462^' ^40) , 

P22,i = 4(54^' +373^' +495^' -6924f -2980) , 

A 1,0 ~ P22,l ' 
Pl2,0 ~ P\2,\ 5 

P22,^ ~ Pn,\ ' 
Ai,-i ~ Pn,\ -> 
A2.-1 = -8(9<?' -142,^' -ISO*?" +1821<^' -670) , 

P22,-1 ~ Pl\,l' 



These equations have the appropriate hmits when 
either £1 or £2 is zero. The mode average is 






(66) 



8(341 -l<+90^: -2434^^+10) 
1125(4^2)^ 

The correctness of Eqs. (64) and (65) was tested by 
calculating the modes of a triaxial ellipsoid using the 
fmite-element method. The parameters £1 and £2 were 
varied, with the ratio held constant at £/£2 = 2. This 
choice corresponds to a uniform splitting of the triplet 
(the case considered in Ref [10]) as shown in Fig. 5. 

Figures 5 and 6 show the close agreement of the 
eigenvalues determined with the finite-element method 
and Eqs. (64) (65). 

5. Concluding Remarks 

The formalism derived in this article can be applied, 
in principle, to arbitrary quasi-spherical resonators 
whose shape can be represented by Eq. (3). Section 3.2 
lists the general principles that determine the possible 
contributions to the general series (32). Once the possi- 
ble terms are identified, the use of symbolic algebra 
software can be used to calculate the terms. For increas- 
ingly complex shapes, this process should be pro- 
grammed so as to minimize human error. 
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Fig, 5. Fractional eigenvalue perturbations for the 11 m modes: the 
points are numerical results determined with the finite-element 
method, the lines represent Eqs. (64)-(65) for an ellipsoid with 
£2 = £il2. The key identifies the lines by the value of the index m. 
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Fig, 6. The differences between the values ofk^^ determined with the 
finite-element method for the 1 1 m modes and the predictions k\^^ of 
Eqs. (64)-(65) for an ellipsoid with £2 = £il2. The lines proportional 
to £l were fit to the plotted points. The key identifies the lines by the 
value of the index m. 
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6. Appendix. Evaluation of Infinite Sums 

The sum S^^f defined by Eq. (25) can be evaluated 
explicitly For the case f ^ 0, consider the contour 
integral 



"L 



1 



^^-S^^-r(r+i)y;(^) 



d^ , (67) 



where C is a rectangular contour with comers at 
(± Xj^, ±yo), with yQ>0 and Xj^ sufficiently large that A'^ 
zeros ofjf(^) he along the positive real axis within 
C. The integrand is bounded on C, so the integral 
approaches zero as N -^ oo-^ the sum of all residues 
within the contour is also zero in this limit. 

The integrand has poles at ±^^„, y^'(/+l), and 

±^^,^ for v= 1,2,... A^. When fi^i the poles are all 
first order, with residues 



n±u- 



1 



jt (L ) 



H^W^))-^rJ^, 



y 



and 



7^(±40=-^ 



htv 



^i-^}.^f^.-i'(i'+i) 



, (68) 



(69) 



(70) 



The sum of all residues 7l(±^f^) is -25^^-, so the condi- 
tion that the sum of all residues is zero implies 



or 



-s,„,+n(^,„)+n(^i'(i'+i))=o, (71) 



*^fit,P' ■ 



(72) 



When ^' = ^ the pole at ^^„ is second order, with residue 



H±L) = - 



4[a-£(i + l)J' 



(73) 



and the sum of the series is 



'^int ~ 






(74) 



Next consider the case £'=0. The second factor in 
the integrand of the contour integral is then unity; the 
residues are 



n±L) = 



1 



and 



y'o (4« ) 

1 



, ^^0, (75) 



^(±4v) 



hit! hfv 



(76) 



There is a single pole at ^qi = 0, so the sum of all 
residues is 



X7^(±^,v)+27^(+4„) 



2 ^ 7o(4) 



which implies 



hln hln Jo ySin ) _ 



(77) 



-+- 



2^'« '^^inhi^ln) 



(78) 



When ^ = the (distinct) poles at ±^o« ^^e second order; 
the residues are -3/(4^ qJ, so the sum is 



1 



4£' 



(79) 
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